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(57) A system and method for color characterization 
and transformation obtain color data representing out- 
put of a color imaging system (56; 58), and convert the 
color data using a color space having a white reference 
vector (60) that is adjusted during the conversion. The 
white reference vector can be adjusted according to 
intensities of the color data being converted. Adjust- 
ment of the white reference vector serves to avoid non- 
uniformities for color imaging systems having different 
imaging bases, and thereby eliminates, or at least 
reduces, the amount of empirical adjustment necessary 
to obtain an acceptable visual match between the color 
imaging systems. 
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Description 

The present invention relates generally to color imaging technology and, more particularly, to techniques for color 
characterization and transformation. 

5 Since the introduction of the CIE (Commission International de TEclairage) color measurement system in the early 
1930's. many different color spaces have been proposed for different applications. A color space, also referred to as a 
color "metric," is essentially a coordinate system by which a color can be quantified. 

A color space can be used to characterize tiie color output of a color imaging system relative to otiier color imaging 
systems. The process of "characterizing" a color imaging system generally refers to tfie calculation of a color response 

10 function for tiie particular color Imaging system using the coordinates of the color sp>ace. The color response function 
is the mathematical correlation between a range of input color values and measured color values ot}tained for output 
generated by the color imaging system in response to such input color values. 

A color space also facilitates transformation of color values to produce a color output in one color imaging system 
that matches tiie color output of a different color imaging system. In terms of imaging science, an "ideal" color space 

15 would allow a user to calculate a color mapping between different color imaging systems and to achieve an acceptable 
color appearance match between such systems without the need for subjective or empirical adjustment by an observer. 
The term "color appearance" refers to tiie psychophysical response of an observer to a color under given levels of illu- 
mination and under given viewing conditions. 

Proposed color spaces differ primarily with respect to the parameters expressed on the coordinate axes and the 

20 manner in which such parameters are calculated. Common to all color spaces proposed for CIE color systems, how- 
ever, are the CIE Standard Observer functions. The CIE Standard Observer functions are based on color matching 
functions and result in a unique set of tristimulus values XYZ for any color measured under specified conditions. The 
tristimulus values XYZ are calculated from the spectral output of eitiier an additive or subtractive color system convo- 
luted with the response function of either a 2 degree or 10 degree Standard Observer. In the case of reflective hard 

25 copy, the spectral reflectance curve is typically convoluted wHh a standard illuminant to estimate the expected spectral 
output of the reflective color. 

One color space that has been the subject of a number of proposals is the CIELAB color space. In this color space. 

L* represents lightness, a* represents redness-greenness, and b* represents yellowness-blueness. According to the 

von Kries chromatic-adaptation transform, the proposed L*a*b* color spaces make use of white reference tristimulus 
30 data. A description of the von Kries chromatic-adaptation transform can be found in Gunter Wyszecki and W.S. Stiles. 

Color Science: Concepts and Metiiods, Quantitative Data and Formulae, section 5.12, John Wiley & Sons, Inc., 1982. 

In general, the von Kries chromatic-adaptation transform involves dividing tiie tristimulus values XYZ ot)tained for a 

color produced by a particular color imaging system by white reference tristimulus values for the system. For example. 

the X, Y. and Z tristimulus values of the color under study can be divided, respectively, by tiie X, Y. and Z tristimulus 
35 values for a perfectly diffuse white reflector. The basic concept of the Von Kries approach is to define all colors, botii 

neuti-al and chromatic, relative to tiie "white reference" representing the XYZ tristimulus values of the perfectiy diffuse 

white reflector. 

Perhaps the most well known of the proposed color spaces is tiie CIE 1976 CIELAB color space. The equations for 
the CIELAB color space are as follows: 

40 

L* = 116(YA'n) -16 [1] 

a* = 500 [(X/X„) - (Y/Y ] [2] 

45 b* = 500 [Qf/y „) - (Z/Z „) ] [3] 

where Xn, Yp. and Z^ are the tristimulus values of a perfectly diffuse white reflector under specified viewing conditions. 
The viewing conditions are determined by (1) the illuminant. e.g.. D50, and (2) the Standard Observer (2* or 10**). 
Another proposed L*a*b* color space ("ICPF LAB"), described by the International Color Consortium in section 2.5 

so of International Color Profile Format (ICPF) Proposal No. ISO/TC 1 30/WG2 N, June 1 0. 1 994. uses X^wYmw^mw "eu 
of XpYpZn. where "mw" denotes the media white reference for a particular imaging system. The media white reference 
refers to either the color of an imaging substrate or the white point produced by a color display This proposed color 
space is based on the ColorSync standard developed by Apple Computer. In this color space, the imaging substrate or 
display serves as tiie white reference instead of a perfectly diffuse white reflector. Witii a color printing system, for 

55 example, the white reference is tiie color of the subtractive printing substi-ate. The media white reference ttierefore rep- 
resents tiie tristimulus XYZ values for tiie unimaged printing subsfrate. With a color display, tiie white reference is ttie 
white point produced by an additive phosphor screen or liquid crystal matrix. 

One color space ("HUNT LAB"), proposed by R.W.G. Hunt, in "Revised colour-appearance model for related and 
unrelated colours," Color Research and Application, volume 16. 1991 . pages 146-165, provides an improved nrKXjel that 
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addresses inconplete chromatic adaptation. The HUNT LAB color space recognizes that chromatic adaptation does 
not completely account for changes due to the white reference of a particular color imaging system. The HUNT LAB 
color space seeks to take such changes into consideration by producing a nrxxjifted set of white reference tristimufus 
values using a complex transformation. 

5 In "Image Color-Appearance Specification Through Extension of CIELAB." Color Research and Application, vol* 

ume 18, number 3. June, 1993, pages 178-190. Mark D. Fairchild and Roy S. Berns propose a variation of the HUNT 
LAB color space, known as "RLAB." The RLAB color space nKXiifies the scaling factors 1/Xn. 1/Yn. I/Z^ according to 
whether the "white reference" is slightly off-white or significantly different from neutral. Spedficatty. the RLAB color 
space introduces a fixed correction factor p calculated as a function of Xp, Y^. and Z^. where X^, Yn. and denote the 

10 "media white" values. In the RLAB color space, the tristimulus values of XYZ are first trarsformed via matrix multiplica- 
tion into "long, middle, and short wavelength" tristimulus values L^ M^, and S^. The equation for p is based on Y^ and 
on Ln, Mn, and Sn, which in turn are functions of Xp, Yp. and Zr. In the RLAB color space, the ratios X/Xn, Y/Y^ and 
Z/Zn become (Li Pl/LJ. (Mi Pw/Mn). (Si Ps/Sp). 

Various comparative tests using hurran observers have been performed to debate the efficacy of the atxTve color 

75 spaces and other color spaces. The existing color spaces discussed above generally are helpful for characterizing rel- 
ative color within a particular color imaging system, given a particular imaging media and a particular set of viewing con- 
ditions. Such color spaces are less helpful, however, in characterization or color transformation between two different 
color imaging systems. One problem not addressed by the existing color spaces is the effect of variation in the white 
references for different color imaging systerr^ due to differences in the imaging substrate or display v/hite point. 

20 As an illustration, the CIELAB color space has been observed to produce nonuniformities when used to match color 
proofing systems having different imaging substrate white references. In particular, the CIELAB color space has been 
observed to produce severe non-uniformities when mapping from a substantially white imaging base to a somewhat 
blue-shifted imaging base. The ICPF color space has been observed to produce uniform mapping in light shades of 
color, but less uniform mapping with more intense colors. Due to such nonuniformity, after a transformation function has 

25 been generated to match L*a*b* values between two color imaging systems, it ordinarily will be necessary for a human 
operator to engage in significant empirical adjustment to obtain an acceptable visual match. The necessary empirical 
adjustment is extremely labor-intensive arvj often can require more than a week for completion. 

In view of the disadvantages associated with existing color spaces, there is a need for an improved color space. In 
particular, there is a need for a system and method for color characterization and transformation that uses an improved 

30 color space capable of avoiding nonuniformities for color imaging systems having different imaging bases. Such a sys- 
tem and method would eliminate, or at least reduce, the amount of empirical adjustment necessary to obtain an accept- 
able visual match between the different color imaging systems. 

In a first embodiment, the present invention provides a system arKi method for characterizing a color imaging sys- 
tem. The system and method, in accordance with this first emt)odiment of the present invention, obtain color data rep- 

35 resenting output of the color imaging system, and convert the color data using a color space having a white reference 
vector that is adjusted during the conversion. 

In a second embodiment, the present invention provides a system and method for characterizing a color imaging 
system. The system and method, in accordance with this second embodiment of the present invention, obtain first color 
data representing output of the color imaging system, convert the first color data to second color data using a color 

40 Space, and adjust a white reference vector of the color space during the conversion according to intensities of the first 
color data. 

In a third embodiment the present invention provides a system and method for performing a color transformation. 
The system and method, in accordance with this third embodiment of the present invention, obtain first color data rep- 
resenting output of a first color imaging system, obtain second color data representing output of a second color imaging 

45 system, convert the first color data to third color data using a color space, convert the second color data to fourth color 
data using the color space, adjust a white reference vector of the color space during the conversion of the first color 
data according to intensities of the first color data, adjust the white reference vector of the color space during the con- 
version of the second color data according to intensities of said secorxl color data, and map the fourth color data to the 
third color data to generate color transformation data. 

so The invention will be described in detail in connection with the drawings, in which: 

Fig. 1 is a functional block diagram of a system for characterizing a color imaging system, in accordance with a first 
emtxxfiment of the present invention; 

FigT'2 is a functional block diagram illustrating the functionality of a system and metfxxi for characterizing a color 
55 imaging system, in accordance with the first embodiment of the present invention: 

Fig. 3 is a functional block diagram of a system for performing a color transtormation. in accordance with a second 
embodiment of the present invention; 

Fig. 4 is a functional block diagram illustrating the functionality of a system and method for performing a color trans- 
formation, in accordance with the second entfxxjiment of the present invention; 
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Hg. 5 is a color space plot comparing relative color responses of a color imaging system and target color imaging 
system after application of a color transformation generated in accordance with the second enijodiment of the 
present invention; 

Rg. 6 is a color space plot comparing relative color responses of the color imaging system and target color imaging 
5 system represented in Fig. 5 after application of a color transformation generated using the CIELAB color space; 

and 

Rg. 7 is a color space plot comparing relative color responses of the color imaging system and target color imaging 
system represented in Fig. 5 after application of a color transformation generated using the ICPF LAB color space. 

10 Rg. 1 is a functional block diagram of a system 1 0 for characterizing a color imaging system, in accordance with a 
first embodiment of the present invention. This characterization, also commonly referred to as "profiling," can be used 
to analyze a single color imaging system, or as a basis to transform the color response of the color imaging system to 
match the color response off another "target" color Imaging system. The system 10 of Fig. 1 may be used to implement 
the steps of a method for characterizing a color imaging system, in accordance with this first embocOment of the present 

15 invention. Accordingly, the method of the first embodiment of the present invention will be described in terms of the 
functionality of system 10. 

As shown in Fig. 1 . system 10 includes a processor 12 that executes a software application program configured to 
carry out the steps of a method for characterizing a color imaging system, in accordance with this first embodiment of 
the present invention. The processor 12 can be realized, for example, by a personal computer such as an Apple Mac- 

20 intosh^" or an IBM PC, or by a computer workstation. Alternatively, processor 12 could be realized by a microprocessor 
that accesses a read-only memory (ROM) into which the application program is loaded. The application program could 
be embedded in a color management software package, such as that provided with the 3M Rainbow^" color proofing 
system, commercially available from Minnesota Mining & Manufacturing Company, of St. Paul, Minnesota. The imple- 
mentation of system 1 0 and the method via software affords flexibility in development and modification. The system and 

25 method alternatively can be implemented by an integrated logic circuit, however, for increased color processing speed. 
In addition to processor 12. system 10 includes means for obtaining color data representing the output of a color 
imaging system. Examples of color imaging systems to which system 10 and the method of the present invention can 
be applied include color printing systems, color display systems, and color projection systems. As shown in Fig. 1. the 
color data can be obtained directly from a color imaging system via a color detector 1 4, or by accessing a color data file 

30 stored in a color data memory 16. With a printing system, for example, color detector 1 4 may be configured to measure 
color values for color patches 1 8 formed on a printing substrate 20. as shown in Fig. 1 . With a color display, color detec- 
tor 14 may be configured to measure color values produced on a phosphor screen or liquid crystal matrix. Similarly, 
color detector 14 may be equipped to measure color values produced by a color projection system. As one illustration, 
color detector 14 could be configured to capture a scene or an animated sequence from a display or projection system. 

35 and generate color values representative of the captured imagery. 

The color detector 14 may comprise, for exannple. a color measurement system such as a Gretag'" SPM 50 color 
measurement device, commercially available from Gretag, Inc.. of Regensdorf, Switzerland, or a densitometer, such as 
an X-Rite color densitometer, commercially available from X-Rite, of Grandville. Michigan. For display or projection sys- 
tem applications, color detector 14 alternatively could comprise a video camera or digital camera. The color data 

40 obtained by color detector 14 may be loaded into color data memory 16 as a color data file, as indicated by line 22. 
Alternatively, the color data obtained by color detector 14 could be directly loaded into a memory associated with proc- 
essor 12, as indicated by line 24. The processor 12 can access the color data file stored In color data memory 16, as 
indicated by line 26, to obtain color data previously detected by color detector 1 4. The color data memory 1 6 may store 
several color data files for a variety of different color imaging systems. Thus, a system user can direct processor 1 2, via 

45 a user interface associated with system 1 0, to select one of the various color data files for purposes of characterizing a 
particular color imaging system of interest to the system user. 

The color data obtained from either color detector 1 4 or color data memory 1 6 preferably represents GIE XYZ tris- 
timulus values for each of a variety of color outputs generated by a color imaging system under study. Alternatively, the 
color data can be converted to CIE XYZ tristimulus values. In addition, the color data preferably represents color out- 

50 puts distritjuted across the color gamut of the color imaging system to thereby provide a broad sampling for purposes 
of color characterization. Thus, in the case off a printing system, color patches 18 represent various combinations and 
gradations of colors formed by dye transfer or ink deposition. The CIE XYZ tristimulus values represent the relative 
amounts of primary color stimuli required to match colors within a CIE color system. The relative values of XYZ are influ- 
enced by the power distribution of the illuminant e.g.. D50. and the CIE Standard Observer function, e.g.. 2° or 10^ As 

55 an afternative to CIE XYZ. the color data can be. for example, in the form of RGB data. CMYK density data, or other 
device dependent data. 

To characterize the color imaging system, in accordance with the present invention, processor 12 converts the color 
data obtained from either color detedor 1 4 or color data memory 1 6 using a color space having a white reference vector 
that is adjusted during the conversion. Specifically, processor 12 converts the color data into second color data and 
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adjusts the white reference vector of the color space during the conversion according to intensities of the original color 
data. In the context of the present invention, intensity refers to degree of difference of a color value from the media white 
or the white point of the color imaging system. The term intensity may refer to both chromatic and achromatic colors. 
By adjusting the white reference vector, processor 12 produces a modified L*a*b* color space. In accordance with the 

5 present invention, the modified L*a*b* color space will be referred to herein as L*a*b* color space. As shown in Fig. 1 , 
processor 12 converts the color data obtained from either color detector 14 or color data memory 16 to produce the 
L*a^^ color data and stores the resulting L*a*b* color data in an L^a^)* memory 28, as indicated by line 30. 

The term "white reference vector" refers to the vector defined by the tristlmulus values X^. Y^. obtained for a 
reference associated with the color imaging system. In the CIELAB color space, the white reference vector is fixed and 

10 comprises tristimulus values X^. Y^. for a perfectly diffuse white reflector. In the ICPF LAB color space described 
by the International Color Consortium in section 2.5 of International Color Profile Formal (ICPF) Proposal No. ISO/TC 
130/WG2 N. June 10, 1994. the white reference vector is fixed and comprises "media white" tristimulus values X^w. 
Yrtw Z(nw tor Imaging base associated with the particular color imaging system. With a color printing system, for 
example, the imaging base is the printing substrate to which colorants will be applied to form an image. With a color 

15 display or color projection system, the imaging base is the white point produced by the display or projection system. 

In the HUNT LAB color space described by R.W.G. Hunt, in "Revised colour-appearance model for related and 
unrelated colours." Color Research and Application, volume 16, 1991, pages 146-165. the white reference vector also 
is fixed. In the RLAB color space, described by Fairchild and Berns. in "Image Color-Appearance Specification Through 
Extension of CIELAB," Color Research and Application, volume 18, number 3, June, 1993, pages 178-190. the white 

20 reference vector is fixed and represented by (Pt/Ln). (Pivi/Mn). (Ps^Sn). where L^, Mp, and are "long, middle, and short 
wavelength" tiistimulus values for media white, and Pl, P|^, and Ps are correction factors calculated according to the 
deviation of the media white from true white. 

Unlike the proposed L*a*b* color spaces discussed above, ttie white reference vector for the L'*'a'*^b^ color space of 
the present invention is not fixed, but rather is adjusted during color conversion. The imaging base of the color imaging 

25 system influences calculation of the white reference, as will be apparent from additional discussion below. However, the 
white reference vector is adjusted during color conversion according to the intensity of the original color data being con- 
verted. Adjustment of the white reference vector based on color intensity, in accordance with the present invention, pro- 
vides more accurate color characterization over the entire range of a color imaging system. In particular, adjustment of 
the white reference vector ensures substantial uniformity of the color characterization for both light color shades and 

30 more intense colors. 

Fig. 2 is a functional block diagram illustrating the functionality of system 10 and a method for characterizing a color 
imaging system, in accordance with the first embodiment of the present invention. The processor 12 can convert the 
color data obtained from either color detector 1 4 or color data memory 16 to second color data using any of tiie existing 
L*a*b* color spaces, such as CIELAB, RLAB. etc.. subject to adjustinent of the white reference vector, in accordance 

35 with the L'^a'^b* color space of the present invention. The CIELAB color space will be used, however, for purposes of 
illustration in describing the functionality of the first embodiment of the present invention. 

With reference to Fig. 2, color data is first obtained for a particular color imaging system by either detecting color 
data produced by the color imaging system, as indicated by block 40. or accessing color data in a color data file, as Indi- 
cated by block 42. The resulting color data, indicated by block 44, comprises an array of CIE XYZ data [(X^, Y^. Z^) . . 

40 . (X|^, Yk. Zk)]. In the array indicated by block 44. ttie vector (Xb, Yb. Zb) represents the color data obtained for the imag- 
ing base associated with the color imaging system. Thus, the vector (Xb. Yb, Zb) represents either the color of an imag- 
ing substrate or the white point produced by a color display or color projection system. The vector (X^. Y^. Z^) 
represents the color data obtained for a maximum color output produced by the color imaging system. More explicitiy. 
the value of X^ is given by tiie value of X which satisfies: 

45 

MAX(ABS{X-Xb)) 

for all values of X measured from all possible colors imaged by the system, where ABS returns the absolute value of 
the difference value in parentheses, and MAX returns the maximum difference value for alt values of X relative to Xb- 

50 Values of Y^ and Z^ can be similarly defined. The array indicated by block 44 also will include a nunrtber of intermediate 
vectors distributed aaoss tiie range of tiie color imaging system. 

In Fig. 2. block 46 indicates the white reference vector (X^. Yp, Z^) for the color imaging system. In the CIELAB.coIor 
space, the white reference vector (Xp. Y„. Zp) ordinarily would represent the color data obtained for a perfectly diffuse 
white reflector. In the L^a^b'^ color space of the present invention, ttie white reference vector (Xn. Yn, Zn) is calculated 

55 as a function of the imaging base vector (Xb. Yb. Zb). and is adjusted according to the particular color data being con- 
verted, as indicated by block 48 of Rg. 2. In particular, tiie white reference vector (Xp. Y^, Zn) is adjusted for each set 
of color data shown in block 44 according to the intensity of the color data. Adjustment of the white reference vector (X^. 
Yn. Zn) results in a new white reference vector (Xn\ Yn'. Zn'). as indicated by block 50. The new white reference vector 
(Xn. Yn . Zn*) tiien is used for conversion of the set of color data. Specifically, as indicated by block 52. the new white 
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reference vector (Xp', Z^') is used as the basis for conversion of the set of color data to modified L*a*b* color data. 
The modified L*a*b* is designated as L'^a'^b"*^ color data in block 54 of Fig. 2. In accordance with the present invention, 
the white reference vector (Xn. Yp. Zn) can be adjusted to form the new white reference vector (Xp . Yn\ Zn*) according 
to the following equations: 

X = X b(1 - SAT(X.Y Z)) + X n * SAT(X,Y Z) [4] 
Yp- = Yb(l - SAT(X,YZ)) + Y„ * SAT(X.YZ) [5] 
Z^. = Zb(1 - SAT(X.YZ)) + Z„ * SAT(X.YZ) [6] 

where: 

SAT(X.YZ) = MAX((X.X„)/(Xk - X„).{Y.Y„)/(Yk - Y„). {Z-Z^)/[Z^ - Z„)). [7] 


where X5. Yt>, Zb are tristimulus values obtained for an imaging base associated with the color imaging system, and 
where X^, Y^, Z^ are tristimulus values of maximum intensity values for XYZ in the color imaging system. The "MAX" 
operator in equation [7] returns the maximum value of the parameters enclosed in parentheses, e.g.. (X-Xn)/(Xk-Xn). 
For color characterization, the adjusted white reference vector (Xp'. Yp'. ZnO can be applied to the CIELAB equa- 
20 tions as follows: 

L* = ^^e(yr/ n")^^ ' ^s [8] 

a* = 500 [(X/X - (YA^ ] [9] 


b* = 500 [(YA^n ) ^'^ ■ i^nl 1 1^01 


Application of equations [8]-[1 0] to the color data obtained for the color imaging system, in combination with the adjust- 
ment of the white reference vector (Xn'. Yn', Zp') as in equations [4]-[7]. produces a set of color data in L+a-^b-" space 

30 that characterizes the color response of the color imaging system. In practice, the values Xk. Y^, 2^ in equations [4]-[7] 
generally can be set to zero for simplicity, because the values of XYZ become very small in the most intense colors, 
e.g.. maximum black. After converting the color data using L+a+b* color space, conventional multi-dimensional interpo- 
lation can be performed to produce a more complete characterization of the color imaging system. 

The application of equations [4]-[7] above to produce the adjusted white reference vector (Xp'. Yp', Zp'), in accord- 

35 ance with the present invention, is based on the following assumptions about the color imaging systems being charac- 
terized. First, it is assumed that the white reference vectors (X^. Yp. Zp) of two different color imaging systems may be 
different even if the illuminant is the same. With color printing systems, for example, the white reference vectors will be 
different for images observed on differently colored printing sit)strates under the same viewing conditions. Similarly, the 
white reference vectors can be different due to different white points produced by different color displays or color pro- 

40 jection systems. Second, it is assumed that the achievable colors of the target color imaging system to be matched are 
generally within the color gamut of the system which is to be transformed in order to nnatch it. even if the white refer- 
ences are not the same. Thus, if the printing substrates or displays on which two images are formed differ by 6 delta E 
according to the 1 976 CIE definition, it can often be the case that the more intense colors can still be matched to 1 delta 
E. 

45 Third, if the first and second assumptions above are correct, then "partial chromatic adaptation" must be reinter- 
preted to mean that von Kries-type scaling, as used in the ICPF LAB color space, is appropriate in the lighter shades 
of color, but is not appropriate in the more intense areas. Thus, regardless of whether one uses CIELAB, ICPF LAB, 
HUNT LAB, RLAB. or any other color space, the white reference vectors can be different for two different imaging sys- 
tems in the less intense colors, i.e., close to the white reference of the imaging base, or "media white." but should con- 

50 verge to be substantially the same in the more intense regions of color space. Even if the media are different, such as 
a display and paper hard copy, substantially the same values for the white reference vectors should be used for the two 
systems for highly intense colors. The L*a*b* color space of the present invention is constructed to cause conversion 
of the white reference vector in more intense regions, as is apparent from equations [4]-[6) above. 

Equations^similar to equations (4]-I7] above can be used for other L*alD* color spaces. For example, the RLAB 

55 color space equations can be modified by adjusting the white reference vector (Lp. Mp. Sp) according to color data inten- 
sity as follows: 

Lp = Lb(1 - SAT(L.M.S)) + L„ * SAT(L.M.S) [111 
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M„ = Mb(1 - SAT(L.M.S)) + * SAT(L.M.S) [12] 

Sn = Sb(1 - SAT(L.M,S)) + S„ * SAT(L.M.S), [13] 

5 where 

SAT(L.M.S) = MAX((L.L„)/(L,-L„). (M-M„)/(M^.M„). (S-S„)/(Sk-S J). [14] 

For simplicity, each of the SATQ equations [7] and [13] above has been represented as a linear function. More complex 
10 functions may be desirable, however, to achieve inproved visual results. For example, the white reference vector can 
be adjusted using the following SAT() function: 

SAT(X.Y Z) -> F(SAT{X,X Z)) [1 5] 

75 where function FQ is an nth order polynomial of the form: 

F(x) = a 0 + a ^x + a2X^+ ... a^x". [16] 

Alternatively, the white reference vector could be ^justed using a SAT() function as follows: 

20 

SAT {XX 2) -> F(L*. C*) [1 7] 

where and are indicative of levels of lightness and colorfulness. respectively, and are based on powers of l/f and 
differences of powers of 1/y. where 3 < 7 < 4.5 depending on surrounding conditions determined according to the HUNT 
25 LAB color space. 

Fig. 3 is a functional block diagram of a system 32 for performing a color transformation, in accordance with a sec- 
ond embodiment of the present invention. The system 32 substantially corresponds to system 1 0 of Fig. 1 . For example, 
system 32 includes a processor 12. a color detector 14. and a color data memory 16. The functionality of system 32 
also substantially corresponds to that of system 10 of Fig. 1. However, system 32 characterizes two or more different 

30 color Imaging systems and calculates a mapping between the color imaging systems based on the characterizations. 
The mapping can be used to generate a color transfomiation table 34, as indicated by line 36. The color transformation 
table 34 then can be used to produce a color response on one color imaging system that visually matches the color 
response on the other, "target" color imaging system. 

In accordance with this second embodiment of the present invention, color detector 14 is configured to obtain color 

35 data representing the output of two or more color imaging systems. As shown in Fig. 3, for example, color detector 14 
may be configured to detect color data from color patches 1 8 formed on a first printing substrate 20 by a first color imag- 
ing system and color patches formed on a second printing substrate 38 by a second color imaging system. With differ- 
ent color monitor displays, color detector 1 4 may be configured to measure color values produced on phosphor screens 
or liquid crystal matrices associated with the displays. Similarly, color detector 14 may be used to measure color values 

40 produced by a different projection systems. The color data for the different color imaging systems alternatively can be 
obtained by accessing color data files stored in color data memory 16. 

As in the first emt)odiment. the color data obtained from either color detector 1 4 or color data memory 1 6 preferat>ly 
refxesents CIE XYZ tristimulus values for each of a variety of color outputs generated by the different color Imaging sys- 
tems. To perform a color transformation, in accordance with the present Invention, processor 12 converts the color data 

45 obtained for each of tiie different color imaging systems using the L*a*b* color space of the present invention, In which 
the white reference vector is adjusted during the conversion. Specifically, processor 1 2 converts first color data obtained 
for a first color imaging system using a color space having a white reference vector that is adjusted according to inten- 
sities of the first color data, and converts second color data obtained for a second color imaging system using a color 
space that is adjusted according to intensities of the second color data. As in the first emljodlment, adjustment of the 

50 white reference vector during conversion of the first color data and the second color data produces modified L*a*b* 
color spaces for the respective data, referred to herein as L'^a^b'^ color spaces. 

Fig. 4 is a functional block diagram illustrating the functionality of system 32 and a method for performing a color 
transformation, in accordance with the second embodiment of the present invention. As shown in Fig. 4, color data is 
obtained for a first color imaging system, as indicated by block 56. and for a second color imaging system, as indicated 

55 by block 58. The resulting color data for the first color imaging system conprises an array of CIE XYZ data [(X^i . Y^-i . 
Zbi) • ■ (Xki. Yki. Zki)], whereas the color data for the second color Imaging system comprises an array of CIE XYZ 
color data [(X^, Y^a. A^) • • • (Xk2. Y^g. 2^)] . As in the first embodiment, the vector (X^i , Y^^ , Zbi) represents tiie color 
data obtained for the imaging base associated with the first color imaging system, and the vector (Xb2. Yb2. Ai2) repre- 
sents the color data obtained for the imaging base associated with the second color imaging system. Similarly, the vec- 
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tor (Xki. Xki, Z^i) represents the color data obtained for a maximum intensity color produced by the first color imaging 
system, and the vector {X^, V^, Z^) represents the color data obtained for a maximum intensity color produced by the 
second color imaging system. Again, the value of X^i is given by the value of which satisfies: 


for all values of measured from all possible colors imaged by the system, where ABS returns the absolute value of 
the difference value in parentheses, and MAX returns the maximum difference value for all values of X^ relative to X^. 
Values of Y^i, Z^i. X|^2, Yk2. and can be similarly defined. 

10 In Fig. 4. block 6o' indicates the white reference vector (Xn. Yp. Zp) for each of the color imaging systems. In the 
L'*'a*b* color space of the present invention, the white reference vector (Xp, Yp. Zp) for each color imaging system is 
calculated as a function of the imaging base vector (Xb. Y^. ZQ for the respective color imaging system, and is adjusted 
according to the particular color data being characterized, as indicated by block 62 of Fig. 4. In particular, the white ref- 
erence vector (Xp. Yp, Zp) for the first color imaging system is adjusted according to the intensity of the color data incD- 

15 cated in block 56. whereas the white reference vector (Xp. Yp, Zp) for the second color imaging system is adjusted 
according to the intensity of the color data indicated in block 58. Adjustment of the white reference vector (Xp, Yp, Zp) 
for each color imaging system results in new white reference vectors (Xpi*. Ypi', Zpi') and (Xn2', Yp2\ Zpa'). indicated by 
blocks 64 and 66, respectively. The new white reference vectors (Xpi*. Ypi*. Zpi') and {Xt^\ Ypa'. ZpaT are used for char- 
acterization of the color data for each color imaging system by conversion using the L'^a^b'^ color space, as indicated 

20 by block 68. The conversion results In new L'^a"^* color data for each of the color imaging systems, as indicated by 
blocks 70 and 72, respectively 

As in the first embodiment, the white reference vector (Xp. Yp. Zp) can be adjusted according to equations [4]-[1 7] 
above during conversion to form the new white reference vectors (Xp^', Ym'. Zn^^ and (Xp2 , Yn2'. Zpa'). After the L*a*b* 
color space data has been produced, multi-dimensional interpolation is performed to complete the characterization of 

25 each of the different color imaging systems. As indicated by block 74, processor 12 calculates a mapping between the 
two different color imaging systems based on the interpolated L^a*b* data. With different color printing systems, for 
example, the mapping may take the form of CMYK -> C'M*Y*K'. The processor 12 can generate a color transformation 
table, as indicated by block 76. defined by the calculated mapping for each of a plurality of color intensity levels. The 
color transformation table can be used to match the output of the first color imaging system to the output of the second 

30 color imaging system. 

The following non-limiting example is provided to further illustrate a system and method for characterizing a color 
imaging system, as well as a system and method for performing a color transformation, in accordance with the first and 
second embodiments of the present invention. 

35 EXAMPLE 

A first 6x6x6x6 grid of color patches was generated by a 3M Rainbow^" desktop color proofing system. A second 
6x6x6x6 grid of color patches was generated using a 3M Matchprint^" III commercial color proofing system. The grids 
were measured by a spectrophotometer and the resulting data stored in a color data file. The grids were ordered as a 

40 matrix of increasing values of CMYK. In the matrix, all permutations of six gray levels were generated. 

Spectral data obtained from each patch was converted into CIE XYZ color space using a 2*^ standard obsen/er and 
D50 illuminant. At this point, the CIELAB color space would dictate the use of equations [1], [2], and [3] above to convert 
the XYZ tristimulus values to L*a*b* color space. The RLAB color space would convert the CIE XYZ values to L^a^b'^ 
using several intermediate matrix multiplications. 

45 Instead, the L+a"*^b"*" color space of the present invention was used, as disclosed herein, by the application of equa- 
tions [4]-[10] to the CIE XYZ color data. Thus, the values of XpYnZp in the CIELAB equations [1]-[3] were replaced by 
Xp'Yp'Zp, as defined in equations [4]-[7]. The new L+a"^b"^ color space was then used to calculate values of L*a*b* for 
the first and second grids generated by the 3M Rainbow^" and 3M Matchprint^" proofing systems, respectively. 

Multi-dimensional interpolation was used to aeate a color profile for each of the two proofing systems, and a CMYK 

so -> C'M'Y'K* mapping was performed with the 3M Matchprint^" system as the target. An example of a suitable mapping 
technique Is described in Stone et al.. "Color Gamut Mapping and the Printing of Digital Color Images,** ACM Transcripts 
on Graphics, volume 7. number 4, October 1988, pages 249-292. After slight manual adjustments were made to the 
resulting color transformation to minimize errors due to measurement noise, uniformity error, etc. a good visual match 
was achieved. 

55 Rg. 5 is a color space plot comparing relative color responses of a 3M Rainbow^" color proofing system and a tar- 
get 3M Matchprint™ color proofing system after application of a color transformation table generated with the L^a^b"*^ 
color space of the present invention to the 3M Rainbow^" color proofing system. In Fig. 5. reference numeral 78 indi- 
cates the color response of the target 3M Matchprint™ system, whereas reference numeral 80 indicates the color 
response of the 3M Rainbow^ system subject to the color transformation. As shown in the plot, the color transformation 


5 


MAX(ABS(Xi - Xb)) 
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provided good visual resufts. The remaining errors that can be seen in Fig. 5 are due to errors in the mapping that could 
be easily fixed with the application o1 more sophisticated software tools. Nevertheless, the color quality was comparable 
to the quality of color transformations previously achieved only with the aid of latX)rious empirical adjustments. 

Fig. 6 is a color space plot comparing relative color responses of a 3M Raintxjw^" color proofing system and a tar- 

5 get 3M Matchprint^" color proofing system after application of a color transformation table generated using the conven- 
tional CIELAB color space to the 3M Rainbow^" color proofing system. Fig. 6 shows the same spectral data as that 
shown in Fig. 5 for 3M Rainbow^" and 3M Matchprint^^ proofing systerr^ plotted using the CIELAB color space. In Rg. 
6. reference numeral 82 indicates the color response of the 3M Matchprint^" system, whereas reference numeral 84 
indicates the color response of the 3M Rainbow^** system. 

10 Fig. 6 reveals significant differences in the color response for the two systems. The differences are due to the ina- 
bility of the CIELAB color space to maintain uniformity across the color range when the colors of the imaging bases 
associated with the two systems differ significantly. For example, Fig. 6 shows a significant shift toward blue in the data 
obtained for the 3M Rainbow™ system. The shift is due to the blue hue of the printing substrate used in the 3M Rain- 
bow^" system. However, the color of the printing substrate must be used for the white reference vector in the CIELAB 

15 color space. Therefore, the mapping procedure must now endeavor to "fix" the blue shift by adding significant amounts 
of yellow based on empirical observations. 

Fig. 7 is a color space plot comparing relative color responses of a 3M Rainbow^" color proofing system and a tar- 
get 3M Matchprint™ color proofing system after application of a color transformation table generated using the ICPF 
LAB color space to the 3M Rainbow™ color proofing system. In Rg. 7, reference numeral 86 represents the color 

20 response of the target 3M Matchprint™ system, whereas reference numeral 88 represents the color response of the 3M 
Rainbow™ system. As shown in Fig. 7, the colors near the origin show similar values between the two systems. How- 
ever, colors extending out from the origin are significantly yellow shifted. The shift is due to a shift in the entire color 
space toward yellow to compensate for referencing all colors to the white reference vector of the imaging base, which 
is blue-shrfled. 

25 Analysis of the average of the errors between the 3M Rainbow™ and 3M Matchprint^" proofing systems with the 
different color spaces described above reveals a systematic difference using the L^a*^"^ color space of the present 
invention, a significant blue shift using CIELAB. and a significant yellow shift using the ICPF color space. 

The average of the errors between the 3M Rainbow™ and 3M Matchprint™ proofing systems using the L^a'*!)'*' color 
space of the present invention were: 

30 


AL* 

Aa* 

Ab* 

-1.15 

-0.52 

0.18 


35 

The average of the errors between the 3M Rainbow™ and 3M Matchprint™ proofing systems using the CIELAB 
color space were: 

40 


AL* 

Aa* 

Ab* 

-1.30 

0.02 

-1.96 


45 

The average of the between the 3M Rainbow™ and 3M Matchprint™ proofing systems using the ICPF color space 
were: 


AL* 

Aa' 

Ab* 

-0.70 

-1.35 

3.52 


55 As indicated above, there is a significant shift toward blue or yellow in the highlight or high intensity regions using 
the CIELAB and ICPF color spaces when comparing two images which appear to be optimally matched in terms of 
overall color. There is a minimal shift, however, when the L^a^^ color space of the present invention is used. A system 
and method employing the L'^a'^'*^ color space therefore is very effective for performing color transformations between 
different color imaging systems, and especially between different color imaging systems having different imaging bases 
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or white points. 
Claims 

1- A method for characterizing a color imaging system, the method comprising the steps of: 
obtaining color data representing output of said color imaging system; and 

converting said color data using a color space having a white reference vector that is adjusted during the con- 
version. 

2. The method of claim 1 , wherein the color data is first color data, the method further conrprising the steps of: 


converting said first color data to second color data using the color space; and 

adjusting the white reference vector of the color space during the conversion according to intensities of the first 
15 color data. 

3. The method of claim 2, wherein said step of converting said first color data includes the step of converting said first 
color data to L^a^b"^ data and further includes the steps of: 

20 converting said first color data to said L*a*b* data according to the following equation: 

L* = 116(YA'n) ^'^-16, 


a* = 500 [(X/Xn)iy3 - (Y/Yn) . and 


wherein X. Y, Z are tristimulus values for said first color data and Xp. Yp. Zn are white reference tristimulus values 
for said color imaging system; and 

adjusting said white reference tristimulus values Xn. Y^. Zp according to intensities of said tristimulus values X, 
Y. Z. thereby adjusting said wvhite reference vector of said color space. 

4. The method of claim 2. wherein said step of converting said first color data includes the step of converting said first 
35 color data to L^a^b*^ data and further includes the steps of: 

converting said first color data to said L*a*b* data according to the following equations: 


L*= 116(YA^n)^^- 16, 

a* = 500 [(X/X . (Y/Y^) and 

b* = 500 [(Y/Yj'^-(Z/Zj'^ ], 

45 wherein X. Y. Z are tristimulus values for said first color data and X^, Y^, Z^ are white reference tristimulus values 
for said color imaging system; and 

adjusting said white reference tristimulus values Xp, Yn. Z^ according to the following equations: 

50 Xn = Xt,(1 - SAT(X,Y.Z)) + X^ * SAT(X,Y.Z). 

Yn = Yb(1 - SAT(X.YZ)) + Yn * SAT(X.YZ). and 

Z„ = Z^(1 - SAT(X,YZ)) + Z„ * SAT(X.YZ). 

wherein: 

SAT(X.YZ) = MAX((X-X J/(X,-X J,(Y-Y„)/(Yk-Y„). (Z-Z J/(Z,-Z J). 
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wherein Z^^ are tristimulus values for an imaging base associated with sakj color imaging system, and 

wherein X^. Y^. are tristimulus values of maximum intensity values for XYZ in said color imaging system. 

5. A system for inrplementing the method of any of clainis 1 -4. 

6. A method for performing a color transformation, the method comprising the steps of: 

obtaining first color data representing output of a first color imaging system; 
obtaining second color data representing output of a second color imaging system; 

converting said first color data to third color data using a color space having a white reference vector that is 
adjusted during the conversion of said first color data; 

converting said second color data to fourth color data using a color space having a white reference vector that 
Is adjusted during the conversion of said second color data; and 

mapping said fourth color data to said third color data to generate color transformation data. 

7. The method of claim 6, wherein the step of converting said first color data to said third color data includes adjusting 
the white reference vector of said color space during the conversion of said first color data according to Intensities 
of said first color data, and wherein the step of converting said second color data to said fourth color data includes 
adjusting said white reference vector of said color space during the conversion of said second color data according 
to intensities of said second color data. 

8. The method of daim 7. wherein said step of converting said first color data includes the step of converting said first 
color data to L*a*b* data, and said step of ccmverting said secorxJ color data includes the step of converting said 
second color data to L*a*b* data, and wherein: 

25 said step of converting said first color data further includes the steps of: 

converting said first color data to said L'a^"^ data according to the following equation: 

L* = 116(Y,/Y„i)''^-16. 

30 

a* = 500 [(Xi,X„,) ''^ - (Y,/Y„i) ''^l, and 

b* = 500[{Y,/Y„,)''^-(Z,/Z^,)''^ ], 

35 wherein X^, Yi are tristimulus values for said first color data, and Xn^. Yp^, Zpi are white reference 

tristimulus values for said first color imaging system, and 

adjusting said white reference tristimulus values Xpi , Yp^ . Zpi according to intensities of said tristimulus values 
X<|, Y^. Z^, thereby adjusting said white reference vector of said color space, and 

40 said step of converting said second color data further Includes the steps of: 

converting said second color data to said L'^a^b* data according to the following equation: 

L* = 116(Y2A^„2) '^-16. 

45 

a* = 500 [{X^/X^) ''^ - (Yg/Ypg) "^1. and 

b* = 500 [(Y2/Y^) - (Z^/Z^) J, 

50 wherein X2, Y2, Z2 are tristimulus values for said second color data, and Xn2, Yn2. Zn2 are white refer- 

ence tristimulus values for said second color imaging system, and 

adjusting said white reference tristimulus values Xn2. Yp2. Zn2 according to intensities of said tristimulus values 
X2. Y2. Z2, thereby adjusting said white reference vector of said color space. 

55 9. The method of dalm 7. wherein said step of converting said first color data includes the step of converting said first 
color data to L'^aV data, and said step of converting said second color data includes the step of converting said 
second color data to L^a^b"^ data, and wherein: 

said step of converting said first color data further Includes the steps of: 
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converting said first color data to said L^a^'b* data according to the following equation: 

L* = 116(Yi/Y„i)^'^-16. 
5 a* = 500 {(X yX„,) - (Y ,/Y ^^]. and 

wherein X^. are tristimulus values for said first color data, and X„i. Y^, are white reference 
10 tristimulus values for said first color imaging system, and 

adjusting said white reference tristimulus values Xpi. Yni, Zni according to the following equations: 

X ni = X (1 - SAT(X 1 .Y , .Z , )) + X „i * SAT(X ^ . Y , .Z ^ ). 

15 Y = Y (1 - SAT(X , .Y , .Z i)) + Y„, * SAT(X 1 .Y 1 .Z , ). 

and 


20 


Z„i =Zbi(1 -SAT(X,.Yi.Zi)) + Z„, *SAT(Xi.Y,.Zi). 
wherein: 

SAT(X , .V ,.2 1 ) = MAX((X , -X )/(X -X ). (Y , -Y )/(Y -Y ). (Z , -Z )/(Z -Z )). 

25 wherein X^i. Yb^. Z^i are tristimulus values for an imaging base associated with said first color imaging 

system, and wherein Xi^v Yi^i , Z^^ are tristimulus values of maximum intensity values for XYZ in said first color 
imaging system, and 


30 


said step of converting said second color data further includes the steps of: 
converting said second color data to said L*a*b^ data according to the following equation: 

L* = 116(Y2/Y^) ^^^-16. 

35 a* = 500 KXg/Xna)''^ -(yz^nz) '^l 

b* = 500 [(Yg/Y ^) ' (Zg/Z ^) ]. 

wherein X2, Y2. Z2 are tristimulus values for said second color data, and Xn2. Yn2, Zn2 are white refer- 
40 ence tristimulus values for said second color imaging system, and 

adjusting said white reference tristimulus values Xn2. Yn2. Zn2 according to the following equations: 

Xn2 = Xb2n - SAT(X2.Y2.Z2)) + X^g * SAT(X g.Y g.Z g). 

^ = Yb2(1 - SAT(X2.Y2.Z2)) + Y„2 * SAT(X2.Y2.Z2). 

and 


50 


Za2=Zb2(1 -SAT(X2.Y2.Z2)) + Z^*SAT(X2.Y2.Z2), 
wherein: 

SATCXg.Yg.Zg) = MAX((X2-X^)/(X^-X^). (Y2-Y ^)/(Yj^-Y (Zg-Z^VCZ^^-Z^)). 

55 wherein Xb2. Yb2. Zb2 are tristimulus values for an imaging base associated with said second color 

imaging system, and wherein Xk2, Yk2. Zk2 are tristimulus values of maximum intensity values for XYZ in said 
second color imaging system. 

10. A system for inrplementing the method of any of claims 6-9. 
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